THE MODELLING AND ANALYSIS OF 
FRACTIONAL-ORDER CONTROL SYSTEMS IN 
FREQUENCY DOMAIN 



I. PETRAS^, IL. DORCAK^, P. O LEARY^, B. M. VINAGRE^, I. PODLUBNY^ 

^ Department of Informatics and Process Control 
BERG Faculty, Technical University of Kosice 
B. Nemcovej 3, 042 00 Kosice, Slovak Republic 
phone: (+42195) 6025172 
e-mail: {petras, dorcak, podlbn}@tuke.sk 

9 

Department of Automation, University of Leoben 
Peter - Tunner strasse 27, A-8700 Leoben, Austria 
phone: (+43) 3842/4029031 
e-mail: oleary@unileoben.ac. at 

Department of Electronic and Electromechanic Engineering 
Industrial Engineering School, University of Extremadura 
Avda. De Elvas, s/n., 06071-Badajoz, Spain 
phone: (+34) 924289600 
e-mail: bvinagre@unex.es 



Abstract 

This paper deals with fractional-order controlled systems and fractional-order controllers in the fre- 
quency domain. The mathematical description by fractional transfer functions and properties of these 
systems are presented. The new ways for modelling of fractional-order systems are illustrated with 
a numerical example and obtained results are discussed in conclusion. 

Key words: fractional-order controller, fractional-order system, fractional calculus, stability. 



1. INTRODUCTION 

Fractional calculus is a generalization of integration and derivation to non-integer order fundamental 
operator a-Df , where a and t are the limits of the operation. The two definitions used for the general 
fractional differintegral are the Griinwald definition and the Riemann-Liouville definition Q]. 

The idea of fractional calculus has been known since the development of the regular calculus, with 
the first reference probably being associated with Leibniz and L'Hospital in 1695. Fractional calculus 
was used for modelling of physical systems, but we can find only few works dealing with the application 
of this mathematical tool in control theory (e.g. ||^, ^, ||, ||]). 

The aim of this paper is to show, how by using the fractional calculus, we can obtain a more general 
structure for the classical PID controller, for controlled systems with memory and hereditary behaviour, 
and how to model fractional-order systems in the frequency domain. On the other hand, we can analyse 
the fractional-order systems and specify the conditions of stability in Bode's and Nyquist's frequency 
response. 



2. FRACTIONAL-ORDER CONTROL CIRCUIT 



We will be studying the control system shown in Fig.|^, where Gdjuj) is the controller transfer function 
and Ga{ji^) is the controlled system transfer function. 
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Figure 1: Feed - back control loop 



2.1 Fractional-order controlled system 

The fractional controlled system will be represented with fractional model with frequency transfer function 
given by the following expression: 

where 13k, otu (fc = 0, 1, 2, . . . , n) are generally real numbers, /?„>...> /3i > Pq, am > ■ ■ ■ > ai > 
and ak,bk (fc = 0, 1, . . . , n) are arbitrary constants. 

Identification methods ||, ^, |l^ for determination of the coefficients Uk , 6fe and ak , (3k {k = 
0, 1, . . . ,rt) were developed, based on minimisation of the difference between the calculated and exper- 
imentally measured values. The optimisation consists in minimising the square norm of the difference 
between the measured frequency response F{uj) and the frequency response of the model Gs{ji^)- The 
quadratic criterion is given by the expression: 



M 

Q^Y. W^{^ra)\F{uJm) " Gsiju^ra)?, (2) 
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where VF(tJm) is a weighting function and M is the number of measured values of frequencies. 
2.2 Fractional order controller 

The fractional PI'^D^ controller will be represented by frequency transfer function given in the following 
expression: 

U{]u:) T, 
— ; r = A + , 

E{ju;) {jujy 

where A and S are arbitrary real numbers {X,S > 0), K is the proportional constant, Ti is the integration 
constant and Td is the derivation constant. 

Taking A = 1 and ^ = 1, we obtain a classic PID controller. If A = and/or Ti = 0, we obtain a PD^ 
controller, etc. All these types of controllers are particular cases of the PI^D^ controller, which is more 
flexible and gives an opportunity to better adjust the dynamical properties of the fractional-order control 
system. 

The PI^D^ controller with complex zeros and poles located anywhere in the left-hand s-plane may 
be rewritten as 

,((jc.)M0^+" + (2g(jc.)^)/c.„-f 1 

where C is a gain, ^ is the dimensionless damping ratio and aj„ is the natural frequency. Normally, we 
choose 0.9 > ^ > 0.7. When ^ = 1, the condition is called critical damping 

The tuning of PI^D^ controller parameters is determined according to the given requirements. These 
requirements are, for example, the stability measure, the accuracy of the regulation process, dynamical 
properties etc. One of the methods being developed is the method (modification of roots locus method) of 
dominant roots ||^, based on the given stability measure and the damping measure of the control circuit. 
Another possible way to obtain the controller parameters is using the tuning formula, based on gain and 
phase margins specifications. 



Gc(jc.) = ^ - ^ + tAx + TS^)'^ (3) 



Gc[3uj)^C j—j^ , (4) 



3. FREQUENCY CHARACTERISTICS 



The graphical interpretation of the frequency transfer function for different vahies of the angular frequency 
uj in the range lo g (0,oo) is called the frequency characteristic. 

The practical meaning of the frequency characteristics is in the determination of the stability of control 
circuits from the behavior of the frequency characteristic. Stability is an asymptotic qualitative criterion 
of the quality of the control circuit and is the primary and necessary condition for correct functioning 
of every control circuit. It is difhcult to evaluate the stability of a fractional-order control circuit by 
examining its characteristic equation either by finding its dominant roots or by algebraic methods. The 
methods given are suitable for integer-order systems. In our case we use Bode's and Nyquist's frequency 
characteristic to investigate the stability . 

For the frequency transfer of the open control circuit Go{juj) with respect to (|l|) and (||), we have 
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It is known from the investigation theory of stability of regulation circuits that the system is stable 
if the roots of the characteristic equation are negative or have negative real parts if they are complex 
conjugate. This means they are located left of the imaginary axis of the complex s-plane of the roots. 
In the case of frequency methods of evaluating the stability we transform the plane of the roots s into 
the complex plane Go{jio) while the transformation is realised according to the transfer function of the 
open circuit. During the transformation all the roots of the characteristic equation are mapped from 
the s-plane into the critical point (— l,iO) in the plane Go{joj). The mapping of the s-plane into plane 
Go{juj) is conformal, that is, the directions and locations of points in the s-plane are preserved also in 
the plane Go{joj). 

In the Nyquist frequency characteristic to find out the stability of the control circuit, it is necessary to 
investigate the behavior of the curve Go{juj) for cu G (0, oo) relative to the critical point (—1; iO). Based on 
the above it follows that also the fractional-order control circuit is stable, if the frequency characteristic 
passes on the right side of the critical point, when going in the direction of increasing values w. If the 
frequency characteristic passes on the left of the critical point, the circuit is non-stable. The passing of 
the frequency characteristic through the critical point means the circuit is on the border of stability. 

4. ILLUSTRATIVE EXAMPLE 



We give in this section an example of modelling the stable dynamical system by using fractional calculus 
in frequency domain. The fractional-order control system consists of the real controlled system with the 
frequency transfer function 

GsO'w) = o.8(jw)2-2 + 0.5(jw)0-9 + 1 

and the fractional-order PD^ controller, designed on the stability measure 5*4 — 2.0 and damping measure 
^ = 0.4, with the frequency transfer function in the form: 



Gcijuj) = 50.0 + 5.3260'tj) 



(7) 
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Figure 2: Nyquist plot (and details) for (^) and 




Figure 3: Bode plots for and (|^) 



In Fig. 2 and Fig. 3 are showed the frequency characteristics of the fractional-order control system 
consists of the fractional-order controlled object (||) and the fractional-order controller (|^). 

5. CONCLUSION 

The above methods make it possible to model and analyse fractional-order control systems in the fre- 
quency domain. The stability of a fractional-order control system can be investigated via the behavior of 
the frequency plot (see Fig. 2 and Fig. 3). These properties of the fractional-order systems can be used in 
the controller parameters design. The results of this and previous works also show that fractional-order 
controllers are robust. This can even lead to qualitatively different dynamical phenomena in control 
circuits. 
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